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1  Introduction 

The  two  applications  addressed  in  this  report  are 

•  spatial  localization  and  scattering  coefficient  computation  from  the  noisy  echoes  of 
chirp  pulses  interrogating  a  multiple  point  target  environment, 

•  numerical  deconvolution  of  Krueger’s  formula  relating  the  reflectivity  kernel  p  of  a 
dielectric  material  to  the  echo  resulting  from  the  interaction  of  an  orthogonally  incident 
chirp  pulse  with  the  material. 

In  parallel  with  these  applications  we  develop  several  topics  in  the  theory  of  the  finite  Zak 
transform  (FZT).  The  FZT  maps  one-dimensional  signals  (vectors)  onto  two-dimensional 
signals  (matrices).  The  Zak  space  representation  provides  a  time-frequency  image  of  the 
one-dimensional  signal.  Signal  decomposition  and  processing  such  as  convolution,  matched 
filtering  and  the  Fourier  transform  can  be  interpreted  geometrically  and  computed  by  matrix 
operations  in  Zak  space. 

A  key  feature  in  the  Zak  space  approach  not  available  in  a  standard  time-domain  or 
frequency-domain  representation  is  that  digital  signal  processing  can  be  viewed  through 
different  time-frequency  windows.  For  applications  addressed  in  this  report  a  specific  win¬ 
dow  Wo  is  taken,  but  windows  can  be  chosen  to  match  classes  of  waveforms,  special  signal 
decomposition  and  processing  and  application  noise  characteristics. 

The  physical  model  in  the  applications  is  the  transmission  of  a  chirp  pulse  and  the  critical 
sampling  of  its  echo  from  a  multiple  point  target  environment  or  a  dielectric  material.  We 
show  that  the  collection  of  shifts  of  a  critically  sampled  chirp  pulse  is  orthogonal  when  viewed 
through  the  window  Wo-  A  negative  result  of  windowing  is  aliasing  which  will  be  discussed 
in  detail,  but  for  the  present  discussion  will  be  ignored. 

Through  W0  a  critically  sampled  echo  of  a  chirp  pulse  can  be  represented  as  an  orthogonal 
linear  expansion  whose  coefficients  in  the  first  application  are  the  scattering  coefficients  of 
the  targets  and  in  the  second  application  are  samples  of  the  reflectivity  kernel.  Matched 
filtering  reduces  to  computing  the  coefficients  of  an  orthogonal  expansion. 

A  closely  related,  equally  important  result  is  that  the  longer  the  time  duration  of  the 
chirp  pulse,  the  finer  its  critical  samples  resolve  targets  and  sample  reflectivity  kernels. 

The  FZT  is  especially  tuned  for  chirp  pulses.  The  main  result  here  is  that  the  support 
of  the  FZT  of  a  critically  sampled  chirp  pulse  and  its  shifts  are  lines  in  Wo-  The  support  of 
a  critically  sampled  echo  is  a  series  of  lines.  Dechirping  results  in  a  series  of  horizontal  lines. 

Under  certain  conditions  the  windowed  Zak  space  representation  of  a  critically  sampled 
echo  of  a  chirp  pulse  can  be  image  processed  to  reduce  noise  and  clutter.  This  processing  can 
stand  alone  or  be  a  preprocessing  step  to  the  windowed  matched  filter  approach  described 
above. 

As  we  will  see  these  conditions  place  constraints  on  the  number  and/or  values  on  the 
lines  in  the  Zak  space  representation  of  the  critically  sampled  echo.  Essential  use  is  also 
made  of  the  vanishing  off  of  these  lines  in  the  noise-free  case. 

The  conditions  under  which  these  image  processing  methods  apply  are  discussed  mainly 
in  Sections  6  and  14.  In  the  first  application  the  relationship  between  the  product  of  the 
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chirp  time  duration  and  chirp  rate  and  the  number  and/or  spacing  of  targets  is  key.  In  the 
second  application  the  relative  sizes  of  the  chirp  time  duration  and  the  material  relaxation 
time  is  the  crucial  factor. 

In  future  work  we  expect  to  extend  these  results  to  multiple  chirp  and  chirp  pulse  train 
interrogation  of  a  multiple  target  scene  consisting  of  point  scatterers,  perhaps  in  motion,  and 
dielectric  scatterers.  Chirp  parameters,  time-duration,  chirp  rate  and  carrier  frequencies  may 
vary  over  the  application. 

A  second  future  goal  is  to  develop  the  FZT  into  a  waveform  design  tool.  Based  on  the 
work  in  this  report,  we  expect  to  design  discrete  waveforms  directly  in  Zak  space  exploiting 
the  relationship  between  the  geometric  structures  supporting  these  waveforms  in  Zak  space 
and  the  effect  of  linear  shifts  and  other  operations  on  these  geometric  waveforms.  The 
specification  of  specially  adapted  windows  in  Zak  space  to  best  make  use  of  this  relationship 
is  a  key  component  to  the  potential  applications  envisioned. 

In  this  report  we  view  the  FZT  as  a  stage  in  the  Cooley-Tukey  (CT)  fast  Fourier  transform 
(FFT)  algorithm.  At  this  stage,  the  one-dimensional  input  data  is  arranged  as  a  two  (or 
higher)-dimensional  data  set  and  the  algorithm  proceeds  by  acting  on  this  two-dimensional 
data  set.  The  advantage  of  viewing  this  computation  in  this  two-dimensional  setting  is  that 
we  can  analyze  the  computation  geometrically.  We  have  used  this  approach  to  construct  an 
orthogonal  eigen  vector  basis  of  the  A-point  FT. 

The  CT  FFT  is  an  example  of  a  divide-and-conquer  algorithm.  This  important  class 
of  algorithms  is  standard  in  many  DSP  computations.  Typically  in  this  approach  a  one¬ 
dimensional  data  set  is  transformed  into  a  multi-dimensional  data  set  and  the  processing 
proceeds  by  operating  on  this  multi-dimensional  data  set.  For  each  of  these  algorithms  we 
can  exploit  geometric  methods. 

The  following  notation  will  be  used  throughout  this  report.  A  vector  x  G  CN  is  written 


x  = 


x0 

X\ 


—  Xn\ 


0<n<N 


and 


L  XN- 1 
XK  =  [^n]o<n<JC  • 


The  inner  product  of  two  vectors  x  and  y  in  C‘  is 


N- 1 

(X:  y)  =  J2  XnVn » 

n=0 

where  *  denotes  complex  conjugation. 

e(/,  0  <  n  <  A,  is  the  vector  in  CN  having  1  in  the  n-th  component  and  0  in  all  other 
components.  The  set 

{«£  :  0  <  n  <  A} 

is  an  orthonormal  basis  of  C^. 


Dn  and  SN  are  related  by 


F(N)SnF(N)~1  =  Dn, 
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and  by 

DrNS%  =  e2^Sr^DrN,  r,  r1  G  Z. 
The  Appoint  time-reversal  matrix  Rm  is  defined  by 


RnY  = 


Since 


we  have 


Vo 

Vn- i 
L  Vi 

F(N)~l  = 


,  y  £  CN. 


F(N)  =  NFiNy'R. 


■N- 


2  Chirps 


A  chirp  pulse  is  any  signal  of  the  form 


xa(t)  = 


e*i-lt2  e2*ipt ,  0  <  t  <  T, 

0,  otherwise, 


t  e 


T  is  the  time  duration,  7  is  the  chirp  rate  and  /?  is  the  carrier  frequency  of  the  chirp  xa. 
Assume  throughout  that 

N  =  7T2 


is  a  positive  integer.  The  methods  developed  in  this  report  apply  to  a  single  transmission 
of  xa.  In  future  work  we  will  extend  these  methods  to  multiple  chirp  transmissions  as  in  a 
SAR  system  and  to  chirp  trains  perhaps  with  varying  time-duration,  chirp  rates  and  carrier 
frequencies. 

7 T  is  usually  taken  as  an  approximation  to  the  bandwidth  of  xa. 

Sampling  xa  at  the  points 


n  n 
~NT=iT' 


n  G  Z, 


we  form  the  critically  sampled  chirp 


xs(n)  =  xa 


em£e2  0  <n  <  N,  neZ 

0,  otherwise, 


where  /  =  pT. 

If  the  chirp  xa  is  transmitted  at  a  point  and  travels  a  distance  r  to  a  reflector  P,  then 
the  echo  at  the  point  is 

2  r 

OiXa(t  t r  ~  • 

C 
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a  is  a  complex  constant  called  the  scattering  coefficient  of  the  reflector  P.  It  includes 
amplitude  attenuation  as  well  as  reflector  reflectivity.  The  echo  from  multiple  reflectors  Pj, 
0  <  j  <  J,  is 

j- 1 

ea(t )  =  <*(j)xa(t  -  trj),  t  E  ft. 

3=0 

In  an  ideal  model  a  multiple  point  target  return  satisfies 


trj  -  Tlj  jy  )  G  Z,  0  ^  j  <  J, 


and  the  critically  sampled  echo  has  the  form 

/  T\  J~l 

es(m)  =  ea  (m—  =  J2  a(j)xs  (m  -  Uj) ,  me  Z. 
v  iV  7  j=o 

The  ideal  model  will  be  assumed  throughout  this  report  except  in  15.4.  However  the  results  in 
this  report  can  be  easily  extended  to  the  following  perturbation  of  the  ideal  model.  Suppose 
R  is  the  distance  between  an  antenna  and  a  point  taken  as  the  origin  in  an  xyz-coordinate 
system  such  that 

2R  T 

—  =  no— ,  n0  >  0,  an  integer, 
c  N 

If  r  is  the  distance  between  the  antenna  and  a  point  in  the  ellipsoid 

x2  +  y2  +  ^-y-  <  R^z, 

N4R 


then 

where 

For 

we  have 

Since 


2  r  T 

tr  =  ~-  n°Ji  +  e> 


T  1 

M  £  N-c 


n . 


0  <  —  T  -tr  <T,  ne  Z 


Xa{^T  -  tr)  =  e«i^\-2Mhra)TteM:LW*T)  e2vri/3(^)Te-2^_ 


7 


n  —  no 
N 


)t|£|< 


n  —  no  1 
N  c’ 


we  can  up  to  a  high  degree  of  accuracy  set 

g-rri-yc2 e-27ri-y(Xzp.^Te  _  ^ 


>  > 
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and  set 


xa  (JjT  -  tr^j  =  f  =  p t. 

Suppose  0  <  e.  Assuming  the  periodicity  condition 

gffiTV g2iri/  _ 


since  0  <  e  <  ^,  we  have 


xa  (n^  -  tT  j  =  e  2irt/3ey(n), 


where 


y(n0  +  m)  = 


me  Z. 


xs(m),  0  <  m  <  N  —  1, 

<  a;s(0),  m  =  N, 

0,  otherwise 

The  perturbation  from  the  ideal  model  introduces  a  shift  which  will  be  more  fully  described 
in  15.4  and  a  carrier  frequency  dependent  phase  factor. 

The  preceding  discussion  can  be  modified  depending  on  the  degree  of  accuracy  required. 
Suppose  A  >  0  is  such  that  the  approximation 


c  =  1,  0  <  n  <  N, 


is  acceptable.  If  r  is  the  distance  between  the  antenna  and  a  point  in  the  ellipsoid 

y2  j, 


x  +y  +  7-5-  <  —FLA, 


then 


where 


LA.  -  M 
N4R  1 


2  r  2  R  ■ 

tr  —  —  — - b  e, 

c  c 


M  < 

~  N  c 


Arguing  as  before  we  have  the  approximation 

T 


where 


xa\  n—-tr)  =  e  2mPiy{n),  n  e  Z, 


(  xs(m),  1  <  m  <  N, 
xs(0),  m  =  N,  me  Z. 


y(n0  +  m)  = 

[  0,  otherwise, 

Allowing  for  zero  scattering  coefficients  we  can  write  es  as 

(R-1)N-1 

es(m)  =  ^  a{n)xs{m  —  n),  me  Z, 

71=0 
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for  a  sufficiently  large  integer  R.  A  coefficient  of  this  expansion  may  be  the  sum  of  the  scat¬ 
tering  coefficients  of  reflectors  equidistant  to  the  transmitter /receiver.  We  call  this  expansion 
for  es  the  discrete  echo  model  for  the  critically  sampled  chirp  xs.  The  critically  sampled  echo 
es  is  the  input  for  processing.  The  problem  is  to  compute  the  scattering  coefficients  a(n), 
0  <  n  <  (R  -  1  )N,  from  the  input  es.  In  Section  7  we  begin  developing  methods  for  solving 
the  problem  in  the  finite  Zak  transform  framework. 

In  Section  7  we  begin  with  an  x  E  CN  and  define  xRN  E  CRN  by 


,  R  —  1  copies  of  0". 


An  echo  of  xRN  is  any  linear  combination 

(R-1)N-1 

e  ~  X/  a{n)SRNxRN ) 

n=0 


xRN  = 


X 

QN 

nN 


where  Srn  is  the  RN  x  RN  shift  matrix.  The  summation  is  taken  so  that  there  is  no  aliasing 
due  to  periodicity  modulo  RN  and  the  cyclic  shifts  are  linear  shifts.  The  theory  developed 
in  this  report  equally  applies  to  multiple  copies  of  xRN . 

The  spatial  resolution  of  xa  is  usually  taken  as  However,  as  we  will  see,  by  orthogo¬ 
nality  we  can  distinguish  two  targets  whose  distances  from  the  an  antenna  differ  by 

T  _  _c_ 

°N  ~  yT 


and  so  we  will  call  ^  the  resolution  of  the  critically  sampled  chirp  xs.  To  achieve  finer 
resolution  we  can  increase  7 T.  In  future  work  we  will  approach  the  problem  of  finer  resolution 
by  over  sampling.  This  method  is  based  on  sampling  xa  at  the  points 


m 

~KN 


T, 


m  E  Z, 


where  AT  is  a  fixed  positive  integer.  The  resulting  oversampled  chirp 


Xovim)  =  Xa 


emK*Ne2mf  0  <  m  <  KN, 
0,  otherwise, 


can  be  decomposed  into  K  sequences  x$,  0  <  k  <  K, 


x^{n)  =  xm{k  +  nK)  =  a;m,(fc)xs(n)e27rt^ ,  n  E  Z, 


where  xs  is  the  critically  sampled  chirp. 
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3  Discrete  Chirps 

In  the  first  part  of  this  report  we  develop  the  algebra  of  discrete  chirps  in  the  framework 
of  the  finite  Zak  transform  (FZT).  The  results  of  this  part  will  be  the  basis  for  studying 
critically  sampled  echoes  beginning  in  Section  8. 

A  discrete  chirp  is  any  sequence  of  the  form 

x(n )  =  e^Ne2"^,  n  £  Z, 

which  is  periodic  modulo  N, 

x(n  +  N)=x(n),  n  £  Z. 

The  periodic  modulo  N  condition  is  equivalent  to 

N  +  2f  £  2Z. 

Throughout  this  report  we  identify  the  linear  space  of  periodic  modulo  N  sequences  with 
CN  by  setting 

=  [y(n)}0<n<N  • 

Sn  is  the  Appoint  cyclic  shift  matrix.  A  cyclic  echo  of  y  £  CN  is  any  linear  combination 

e  =  E  a(r)sNy- 

r= 0 

If  x  is  a  discrete  chirp,  then  the  cyclic  echo  expansion 

e  =  a{r)SrNx 

T— 0 

is  an  orthogonal  expansion  and  its  coefficients  can  be  computed  from  N  inner  products,  a 
special  case  of  match  filtering. 

If  x  is  a  discrete  chirp,  then  the  multiplication  operator  Mx»  is  called  x-dechirping.  The 
x-dechirping  of  a  cyclic  echo  e  of  x  is 

JV-l 

Mx*e(n)  =  Oi{r)x(—r)e~2'niBE ,  n  £  Z. 

r= 0 

Dechirping  provides  a  second  way  to  compute  the  coefficients  of  a  cyclic  echo  e  of  a 
discrete  chirp  x.  First  x-dechirp  the  cyclic  echo,  compute  its  A-point  FT  and  diagonal 
matrix  multiply  the  result  by  the  inverse  of 

ND  ([x(-r)]0<r<N)  . 


y{  0) 
1/(1) 
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4  Finite  Zak  Transform  (FZT) 

In  this  section  we  introduce  the  finite  Zak  transform  (FZT).  For  each  positive  divisor  L  of 
N  we  define  a  linear  isomorphism 


ZL:CN  —^CL  xCK,  N  =  LK, 

of  CN  onto  CL  x  CK .  If  y  €  C^,  then  ZL y  is  a  two-dimensional  L  x  K  time-frequency 
image  of  y. 

For  a  positive  divisor  L  of  IV  and  y  €  CN  construct  the  vectors  zk,  0  <  k  <  K,  in  CL, 


z  k 


Vk 

Vk  +  K 


=  [yfc+i/c]o<;<L  • 


L  Vk+(L-1)K  J 

zk  is  the  subvector  in  CL  formed  by  the  components  of  y  over  the  index  set 


{k  +  lK  :0  <  l  <L}. 


The  sequence 

z0,  Zi,  . . . ,  Z.K-1 

describes  the  evolution  of  these  subvectors.  If  L  =  1,  we  have  the  standard  ordering  of 
the  components  of  y.  For  L  =  2,  we  have  a  description  of  the  evolution  from  the  even 
components  to  the  odd  components  of  y. 

Identify  CL  x  CK  with  the  space  of  L  x  K  complex  matrices.  Define  the  mapping 
Ml  :  CN  — ►  CL  x  CK  by 


Ml y  =  [z0  zj  •  •  •  z*_i] ,  yeC", 
and  the  mapping  Zl  '■  CN  — >  CL  x  CK  by 

Zl y  =  F(L)Ml y,  y  e  CN. 

ZlY  provides  a  two-dimensional  Lx  K  image  of  the  evolution  of  the  FT  of  the  sequence  of 
subvectors. 

ZL y,  y  £  CN,  is  called  the  L  x  K  FZT  of  y  or  more  suggestively,  the  L  x  K  Zak  space 
representation  of  y.  Different  factorizations  of  N  lead  to  different  Zak  space  representations 
each  depicting  different  time-frequency  aspects  of  y  in  the  sense  discussed.  For  example  Z iy 
is  the  transpose  of  y  while  Z^y  is  the  N- point  FT  of  y.  Throughout  this  report  the  theory  is 
developed  for  a  single  factorization  of  N.  A  future  effort  will  be  directed  toward  combining 
two  or  more  Zak  space  representations  for  the  purpose  of  echo  analysis  and  noise  reduction. 

Zl  is  a  linear  isomorphism  of  CN  onto  CL  x  CK  with  inverse 

Zl1  =  M(L)-1F(L)-1. 
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Up  to  scale  factor  Zl  is  an  isometry. 

|Ziy|2  =  £jy|2,  y  6  C". 

An  orthogonal  basis  in  C'v  is  mapped  onto  an  orthogonal  basis  in  CL  x  CK . 

In  this  section  we  derive  general  formulas  relating  the  FZT  of  a  signal  with  that  of  its 
cyclic  shifts.  Suppose  x  is  a  vector  in  C*v  and  set 

ZLx=[X0X1  •••  XK-!],  XkeCL. 

The  FZT  of  the  cyclic  shift  of  a  signal  is  not  simply  the  cyclic  shift  of  the  columns  of  the 
FZT  of  the  signal.  The  cyclic  shift  has  its  0-th  column  multiplied  by  Dl  and  more  generally 

'  (Sjf  x)  =  DnLZLx,  0  <n<L. 

In  this  case  the  columns  of  ZLx  are  multiplied  by  D£. 

For  x  and  y  in  C^,  the  A-point  cyclic  convolution  is  defined  by 

N—l 

U  =  X  *  y  =  xrSrN y. 

r=0 

5  FZT  and  CT  Factorization 

There  are  several  ways  to  motivate  the  definition  of  the  FZT.  We  choose  the  approach  to  the 
FZT  most  directly  related  to  computation.  The  FZT  can  be  viewed  as  an  implementation  of 
two  stages  of  the  Cooley-Tukey  (CT)  FFT  algorithm.  We  approach  the  CT  FFT  algorithm 
through  factorizations  of  the  A-point  FT  matrix  F(N). 

’  We  continue  to  use  the  notation  in  Section  4.  The  A  x  A  stride  permutation  matrix 
P(A,  K )  is  defined  by 

P(N,  K)y  =  h]0£l<K ,  y  E  C". 

Ik  ®  F(L )  is  the  A  x  A  block  diagonal  matrix 

'  F(L)  0 

IK  <8>  F(L)  = 

0  F(L) 

Set  zi  equal  to  the  A  x  A  matrix 

zl  =  (Ik®F(L))P(N,K). 

In  [7],  we  show 

F(N)  =  Kz-kx  ((/L  <g>  Rk)  P{N,  L)Tn{L ))  zLi 
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where  Rk  is  if-point  time-reversal  and  T^{L)  is  a  diagonal  matrix  whose  diagonal  entries 
are  IV-th  roots  of  unity. 

It  is  convenient  and  standard  to  view  the  iV-point  CT  FFT  as  given  by  this  CT  factoriza¬ 
tion  as  a  multi-stage  process  in  which  a  vector  y  e  CN  is  first  mapped  onto  a  two-dimensional 
array,  operations  are  performed  on  the  two-dimensional  array  and  then  the  two-dimensional 
array  is  mapped  back  to  the  vector  F(N)y. 

CT  FFT  computation  of  F(N) y,  y  G 

•  Form  the  L  x  K  array  Ml y. 

•  Compute  the  L-point  FT  of  every  column  of  MLy. 

ZLy  =  F(L)Ml  y. 


These  two  stages  correspond  to  the  matrix  product  zlY- 

•  Multiply  the  coefficients  of  the  Lx  K  array  ZL  y  by  AT-th  roots  of  unity  as  determined 
by  Tn(L). 

•  Transpose  the  resulting  array  to  form  a  K  x  L  array. 

•  Compute  the  LT-point  time-reversal  of  each  column  of  the  transposed  K  x  L  array. 
The  preceding  two  stages  correspond  to  operating  by  the  permutation  matrix  RnP(N,L). 

•  Form  F(N)y  by  operating  on  the  K  x  L  array  by  KZ^1 . 

The  final  stage  correspond  to  matrix  multiplication  by  Kz^1. 

6  FZT  of  Dechirped  Echo 

Throughout  set  w  =  and  v  —  e2rrti . 

In  this  section  we  compute  the  FZT  of  a  dechirped  cyclic  echo  of  a  discrete  chirp  x 

e  =  E  a(r)SN x 

r— 0 

and  compare  the  result  with  the  FT  of  e. 

Define  the  L  x  K  matrix  F^(r),  re  Z,  by 

Fk(r)  =  [e^  -  ■  •  erL]  . 

We  can  view  F^(r)  as  a  horizontal  line  in  L  x  K  Zak  space  at  position  r  mod  L. 
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Set  y  =  Mx*SrNx,  0  <  r  <  N.  To  compute  the  k-th  column  of  ZLy  form  the  A:-th  column 
of  Ml y 

Moskl  =  x(-r>~rt  [.-*]w  =  xf-r^di., 
and  then  take  its  L-point  FT 

Lx(-r)w~rke$. 

Figure  1  displays  the  Zak  transform  of  dechirped  discrete  chirp  and  its  cyclic  shifts. 
The  two-dimensional  intensity  data  are  displayed  as  log-scaled  images.  The  discrete  chirp 
parameters  are  7  =  .25,  T  —  80,  /  =  40.  N  =  7 T2  =  1600.  The  FZT  parameters  are  L  =  50 
and  K  =  32. 


Figure  1:  Zak  transform  of  dechirped  discrete  chirp 


real  part  imaginary  part 


real  part  imaginary  part 


cyclic  shift  by  10 


We  have  that 

{SrNx  :  0  <  r  <  N} 

is  an  orthonormal  basis  of  CN  and  the  set 

{ZL  (Mx.S^x)  :  0  <  r  <  N} 
is  an  orthogonal  basis  of  CL  x  CK. 

Since  the  support  of  the  L  x  K  Zak  space  representation  of  Mx.  SyyX  is  the  horizontal 
line  at  position  r  mod  L,  if 

r  =  r'  mod  L,  0  <  r,  r'  <  N, 
we  cannot  distinguish  the  vectors 

Mx.5^-x  and  M'x.S'^x 

by  the  supports  of  their  L  x  K  Zak  space  representations.  However  if 


r  ^  r'  mod  L,  0  <r,r'<  N, 
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then  the  supports  of  the  L  x  K  Zak  space  representations  of  these  vectors  are  disjoint.  In 
particular  the  supports  of 

{ZL  (Mx.SrNx)  :  0  <  r  <  L} 

are  pairwise  disjoint. 

Generally  the  support  of  the  L  x  K  Zak  space  representation  of  Mx*  SrNx  determines 
r  mod  L  but  not  r  itself.  To  compute  r  we  must  look  at  the  values  of  Zl  (Mx»5'^x)  on  the 
horizontal  line  at  position  r  mod  L  with  values  given  by  the  vector  v6  CK 

v  =  i*(-r)[m— ' =]w.  (1) 


Write 

r  =  a  +  bL,  0  <  a  <  L,  0  <b  <  K. 
a  is  known  and  b  is  to  be  determined.  Then 

“  —bk 


v  =  Lx(-(a  +  bL))Dj,a(K)  [u 
To  determine  b  we  compute 


0<k<K  ’ 

F(K)D%(K)v  =  cef , 


u  =  e2™*. 


where  c  is  the  complex  constant 


c  =  Nx(—(a  +  bL)). 


b  is  determined  as  the  index  of  nonzero  component  of  the  computation. 

Figures  2  and  3  displays  the  Zak  transform  of  dechirped,  cyclically  shifted  discrete  chirps. 
The  shift  amounts  are  of  the  form  a +  bL,  a  =  11,  and  all  congruent  modulo  L  =  50.  The 
second  row  of  the  figure  are  plots  of  the  values  on  the  horizontal  line  through  a.  The  third 
row  of  the  figure  are  plots  of  the  vector  F{K)Dn{KY  =  ceff.  Note  that  the  position  of 
nonzero  value  corresponds  to  b. 

For  a  cyclic  echo  e  of  a  discrete  chirp  x 

e  =  J2  a(r)SN*> 

r=0 

we  have 

ZL(Mx.e)  =  L  £  a(r)x(-r)Fk(r)DJ(K). 

r=0 

Since  the  set 

{x(— r )Fk(r)DJ,r(K)  :  0  <  r  <  N] 

is  an  orthogonal  basis  of  CL  x  CK,  we  can  compute  the  coefficients  a(r),  0  <  r  <  N,  by 
inner  products. 

There  are  two  important  special  cases  for  applications.  Suppose 

e  =  E  ^rjS^x. 

r= 0 
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Figure  2:  Zak  transform 
cyclic  shift  by  11+2L 


of  dechirped  discrete  chirp 

cyclic  shift  by  11+7L 


We  say  that  e  is  a  cyclic  echo  from  closely  spaced  targets.  The  values  of  Zi  (Mx»e)  on  the 
horizontal  line  at  position  0  <  r  <  L  are  given  by  the  vector  vr  G  CK 

vr  =  a(r)x(—r)  [u>_rA:]  , 

r  w  v  '  L  10<k<K  ’ 

where  we  have  divided  out  L.  For  simplicity  in  the  discussion  assume  the  coefficients  of  e 
are  positive  real  numbers  or  zero.  Then 

a(r ) 
a(r) 

In  the  presence  of  noise  we  may  have  errors  attached  to  the  coefficients  of  |vr|.  We  can 
estimate  a(r )  by  best  estimation  techniques  on  the  K  components  of  the  noisy  |vr|.  Figures 
4-6  illustrates  results  of  this  denoising  method.  Figure  4  displays  the  noise-free  cyclic  echo 
which  is  the  sum  of  5  cyclically  shifted  chirps  with  non-negative  coefficients.  The  shifts  are 
9,  17,  26,  30,  32  and  the  respective  coefficients  are  0.5992,  3.0385,  0.9145,  1.8699,  1.7293. 
Figure  5  displays  the  noisy  echo  obtained  by  adding  noise  to  the  cyclic  echo  in  Figure  4. 

The  Zak  transform  of  the  dechirped  noisy  echo  is  used  as  input  to  image  processing 
algorithm  for  detecting  and  locating  the  horizontal  lines.  Location  is  exact  in  this  case. 
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Figure  4:  Noise-free  cyclic  echo 


real  part  imaginary  part 


Figure  5:  Noisy  echo 


real  part  imaginary  part 


illustrate  estimation  of  arbitrary  scattering  coefficients.  The  discrete  chirp  parameters  are 
7  =  .25,  / 0  =  40 ,  N  =  1600.  Figure  12  displays  the  noise  free  echo  which  is  the  sum  of  4 
cyclically  shifted  discrete  chirps.  The  Zak  transform  parameters  are  L  —  40,  K  =  40. 

The  estimated  shifts  are  accurate  at  10,  16,  23,  29.  In  Figure  16  plots  of  the  absolute 
values  of  the  ratio  of  the  detected  horizontal  lines  and  the  corresponding  vectors  given  in 
(!)• 

Table  3  compares  the  test  values  with  estimated  values  of  the  scattering  coefficients. 

The  FT  approach  essentially  sums  these  errors  to  a  single  noisy  component  by  a  iF-point 
FT  computation. 

The  same  argument  holds  for  a  cyclic  echo  of  a  discrete  chirp  x  of  the  form 

L—l 

e  =  a(r)S^r°x,  r0  some  integer. 

r=0 

Observe  that  the  Zak  space  representation  of  e  does  not  directly  determine  the  first  target, 
ro  must  be  known  either  from  the  application  or  computed  by  the  above  discussion.  In  future 
work  we  will  study  the  use  of  the  information  gained  from  several  size  FZT  computations  to 
remove  ambiguities  of  this  kind. 
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Figure  6:  dB  plot  of  signal  to  noise  ratio 


Figure  7:  Zak  transform  of  the  dechirped  noisy  echo 


Assume  that  the  target  structure  is  sparse  relative  to  the  sampling  in  the  sense  that  there 
are  relatively  few  horizontal  lines  in  the  Zak  space  image  of  the  dechirped  echo  as  compared 
with  L.  Since  we  know  that  the  FZT  of  a  noise  free  dechirped  echo  must  be  supported 
on  horizontal  lines  we  can  search  Zak  spaces  for  the  existence  and  positions  of  these  lines, 
determine  noise  statistics  under  the  knowledge  that  the  Zak  space  representation  of  the 
deichirped  cyclic  echo  must  vanish  off  of  the  lines  and  use  the  noise  statistics  to  correct  for 
noise  on  the  lines  themselves.  In  the  first  line  determination  step  if  there  are  any  variations 
in  the  line  structures,  we  can  correct  these  variations  by  best  line  approximation  methods. 

For  sparse  target  structures  the  FZT  approach  restricts  scattering  coefficient  computation 
to  nonvanishing  lines  in  Zak  space  while  the  FT  approach  requires  computation  for  every 
horizontal  line  in  Zak  space. 

7  FZT  of  Discrete  Chirps 

In  this  section  we  describe  the  FZT  of  discrete  chirps.  Dechirping  is  a  powerful  tool  for  the 
applications  considered  in  Section  6  using  either  the  FT  or  FZT.  However  one  of  the  goals 
of  this  report  is  to  show  how  the  framework  of  FZT  can  be  used  to  analyze  echoes  in  a  more 
general  setting  than  that  obtained  in  Section  6.  For  example  the  framework  should  allow  for 
multiple  chirp  sets  and  chirp  pulses,  perhaps  at  varying  chirp  rates  and  should  include  tools 
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Figure  8:  Noise-free  cyclic  echo 


real  part  imaginary  part 


Figure  9:  Noisy  echo 


real  part  imaginary  part 


for  realizing  linear  shifts  as  well  as  cyclic  shifts.  The  results  of  this  section  provide  the  basis 
for  these  extensions  which  will  be  addressed  in  future  work.  At  the  same  time  these  results 
will  suggest  the  use  of  the  Zak  space  framework  in  waveform  design,  another  topic  in  future 
work. 

Define  the  L  x  K  matrix  Ej<( 0)  by 

Ek( 0)  =  [eo  e£x  •  •  •  e^_1}]  . 

As  usual  is  defined  for  m  mod  L.  We  can  view  E^( 0)  as  the  line  of  slope  1  through  the 
origin. 

Set 

Co  =  F{L)Diag  [x(IK)\0<1<lF(L)~\ 

The  support  of 

Co1Zlx  =  LE^(0)D(xk) 

is  £^(0),  the  line  of  slope  1  through  the  origin. 

Define  the  L  x  K  matrix  Ej^(r)  by 

Ek{r)  =  [er  er-l  ■  ■  ■  er-(K-l)\  »  r  G  Z. 
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Figure  10:  dB  plot  of  signal  to  noise  ratio 


Figure  11:  Zak  transform  of  the  dechirped  noisy  echo 


6  10  is  ao  as  30 


We  can  view  E^(r)  as  the  line  of  slope  1  through  the  point  (r,  0)  in  Lx  K  Zak  space.  As 
with  the  horizontal  lines  F^(r)  in  Section  6,  we  have  that 

ELK(r)  =  Ej}  (r') ,  r,  r'  €  Z, 


if  and  only  if 
and  that 
are  disjoint  whenever 
In  particular  the  lines  of  slope  1 


r  =  r'  mod  L 

Ek(r)  and  Ek  (r>)  i  r,  r'  €  Z, 
r  ^r'  mod  L. 


[ELK{r):  0  <  r  <  L}  , 

are  pairwise  disjoint  in  CL  x  CK. 

As  before,  the  set 

(SrNx) ,  0  <r<N, 

is  an  orthogonal  basis  in  CL  x  CK.  Since  Co  is  a  unitary  matrix  we  have  that  the  set 


{ColZL  (SrNx)  :  0  <  r  <  N} 
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Table  2:  comparison  of  estimated  scattering  coefficients 


test  coefficients 

1.0824 

1.2694 

0.6954 

1.5013 

2.1946 

1.6165 

1.8878 

3.9520 

average  coefficients 

1.0263 

1.2099 

0.7001 

1.5277 

2.1236 

1.5740 

1.8178 

3.9048 

median  coefficients 

1.0968 

1.2717 

0.6831 

1.5904 

2.1976 

1.5031 

1.8072 

3.8628 

Figure  12:  Noise-free  cyclic  echo 


real  part  imaginary  part 


is  also  an  orthogonal  basis  in  CL  x  CK. 

Since  the  support  of 

C^ZL  (SrNx)  =  LELK{r)D  (xK)  DN(K)~r 

is  Eft(r),  the  line  of  slope  1  through  the  point  (r,  0),  the  L  x  K  Zak  space  representations 

{Cq'Zl  (S^x)  :  0  <r<L] 

are  pairwise  disjoint. 

The  discussion  in  Section  6  on  the  use  of  the  FZT  operating  on  a  dechirped  cyclic  echo 
can  be  directly  carried  over  to  the  use  of  the  FZT  operating  on  a  cyclic  echo.  The  horizontal 
lines  Fx(r)  are  replaced  by  the  lines  E^(r)  of  slope  1. 

The  following  condition  is  involved  in  theory  and  numerical  experiments  for  notational 
simplicity.  It  removes  the  necessity  of  multiplying  by  Cq 1  to  get  a  simple  geometric  inter¬ 
pretation  of  results. 

We  say  that  a  discrete  chirp  x  satisfies  the  Lx  K  Zak  space  condition  if 

x(lK)  =  1,  0  <  l  <  L. 

The  L  x  K  Zak  space  condition  is  equivalent  to  the  condition 

l2M  +  2fjE2Z,  0  <1<L,  M  =^r. 

L/  L 

The  FT  and  the  FZT  of  a  dechirped  cyclic  shift  of  a  discrete  chirp  x  result  in  sparse 
outputs,  a  single  nonzero  component  for  the  FT  computation  and  a  single  horizontal  line  of 
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Figure  13:  Noisy  echo 


real  part  imaginary  part 


Figure  14:  dB  plot  of  signal  to  noise  ratio 


values  for  the  FZT  computation.  We  have  just  seen  that  the  FZT  of  x  is  also  sparse  having 
support  a  single  line  of  slope  1.  The  analogy  collapses  when  computing  the  FT  of  x. 

Set 

v  =  F(N)x. 

A  direct  computation  shows  that 

jv-i  2 

va  =  ^2  w **  w^a+^n,  0  <  a  <  N. 

n= 0 

Suppose  for  example  N  —  4  and  /  is  even.  Then 

vq  =  2  m-f, 

ui  =  2, 

V2  = 
v3  =  2. 
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Figure  15:  Zak  transform  of  the  dechirped  noisy  echo 


Figure  16:  Magnitudes  of  the  values  on  the  detected  lines/ vectors  given  by  (1) 


shift  by  10  shift  by  16  shift  by  23  shift  by  29 


8  Zero-padding 

In  preceding  sections,  we  developed  the  general  theory  of  the  FZT  for  vectors  in  CAr  or 
equivalently  for  periodic  modulo  N  sequences  and  applied  the  theory  to  discrete  chirps  and 
their  cyclic  shifts.  The  general  theory  will  be  extended  to  zero-padded  vectors  and  applied 
to  zero-padded  discrete  chirps  and  their  linear  shifts. 

Suppose  x  G  CN  and  R  is  a  positive  integer.  Define  the  zero-padding  xRN  of  x  to  CRN 

by 

X 

qN 

xRN  =  eR  tg>x=  .  ,  R  —  1  copies  of  0A  adjoined  to  x. 

As  discussed  in  Section  2,  the  cyclic  shifts 

S%nxrn,  0  <m<(R-  1  )N, 

are  linear  shifts. 

We  begin  by  describing  the  general  theory  of  the  FZT  for  zero-padded  vectors.  By  making 
use  of  the  special  form  of  zero-padded  vectors  and  the  finer  details  of  the  RN- point  CT  FFT 
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Table  3:  comparison  of  estimated  scattering  coefficients 


test  coefficients 
average  coefficients 
median  coefficients 


1.5740+3.99521 
1. 7173+4. 0692i 
4.3647+4.07591 


7.2175+4. 7044i 
7.3979+4. 6806i 
10. 3068+4. 7747i 


-1. 7549-4. 9605i 
-1. 5806-4. 9725i 
0. 9659-5. 0327i 


3. 1162+1. 0602i 
3.3206+1.16801 
6. 3260+1. 4081i 


algorithm,  we  show  that  the  RL  x  K  FZT  of  xN  consists  of  the  L  x  K  FZT  of  x  over  a 
region  of  RL  x  K  Zak  space  and  filtered  versions  of  this  FZT  over  other  regions. 

The  cyclic  shifts  of  xRN  are  more  difficult  to  handle  in  the  RL  x  K  FZT  framework.  The 
reason  for  this  is  that  the  Fourier  transform  does  not  commute  with  the  shift  operations. 
The  general  theory  of  linear  shifts  is  studied  in  Section  10. 

To  compute  the  RL  x  K  Zak  space  representation  of  xRN  we  will  first  derive  a  formula 
for  computing  the  J?L-point  FT  of  a  zero-padded  vector  uRL,  u  £  CL.  The  main  tool  needed 
is  the  CT  FFT  factorization  for  the  RL-point  Fourier  transform.  Recall  the  following  tensor 
product  notation. 

P(RL,  L)  is  the  RL  x  RL  stride  by  L  permutation  matrix  defined  as  follows.  Suppose 
v  £  CRL.  Form  the  vector  wq  £  C72  by  striding  through  v  with  stride  L. 


w0 


Off-setting  by  1,  form  the  vector  Wi 


vo 

vL 

V(R-l)L  . 


Wi  = 


Vl+ 1 


Continuing  in  this  way,  form  the  L  vectors  in  CR 


w0,  Wi,  ...,  WL_j. 


Define 


P{RL,  L)v  = 


w0 

W] 


L  WL_  j 

Define  the  L  x  L  diagonal  matrix  Drl(L)  by 


t  —  e 


2-ni-. 
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and  the  RL  x  RL  diagonal  matrix  Trl  by 


Tri(L)  =  R^D'rl(L). 

r= 0 


Set 

and 


Crl(L)  =  F(L)Drl(L)F(L) 


h 

Crl(L) 


-l 


c  rl(L)  = 

[  Crl1^)  J 

We  will  use  the  following  CT  factorization  for  F(RL). 

F(RL)  =  P(RL ,  L)  ( IR  ®  F(L))  Trl(L)  (F(R)  ®  IL ) . 
For  the  rest  of  this  report  set 


P  =  P(RL,R)  =  P(RL,L)-1 

and 

C  =  CrL(L)  and  C  =  Crl(L). 


PF(RL)url  =  CF(L)u  = 


F(L)  u 
CF(L)  u 


Cr~'F(L)vl 


Segmenting  PF(RL) uRL  into  R  contiguous  segments  of  size  L,  the  first  segment  is  F(L) u, 
while  the  remaining  segments  are  filtered  versions  of  F(L)  u. 

Segmenting 


PZrlxrn  -  C  ZLx  = 


ZLx 

CZLx 


CR~lZLx 


into  R  contiguous  segments  of  size  L  x  K,  the  first  segment  is  ZlX,  while  the  remaining 
segments  are  filtered  versions  of  Z^x. 


9  Zero-padded  Discrete  Chirps 

We  continue  using  the  notation  P  —  P(RL,  R),  C  =  Crl{L )  and  C  =  C rl(L). 
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The  problem  of  computing  the  RL  x  K  FZT  of  a  zero-padded  vector  xRN ,  x  E  CN,  can 
be  reduced  to  that  of  computing  ZLx  and  then  placing  ZLx  into  the  matrix  construction 

PZrlxrn  =  CZLx. 


We  can  now  compute  the  RL  x  K  FZT  of  the  zero-padded  discrete  chirp. 

Suppose  x  is  a  discrete  chirp.  The  vector  xRN  is  called  a  zero-padded,  discrete  chirp. 
Since  circulant  matrices  commute  we  can  write 


a, 

c0c 


C0CR -1 


In  the  segment  of  RL  x  K  Zak  space  formed  by  the  first  L  rows,  the  support  of  PZRlxrn 
is  ER( 0),  while  in  the  remaining  segments  we  ’fill’  the  segments  through  the  actions  of  Cr, 
1  <  r  <  R. 

Similar  arguments  can  be  used  to  derive  flexible  size  FZT’s. 


10  Windowed  Linear  Shifts 

One  of  the  advantages  of  representing  an  echo  in  Zak  space  is  that  the  echo  can  be  viewed 
through  different  time-frequency  windows  in  Zak  space.  Knowledge  of  the  shifts  making  up 
the  echo  can  be  used  to  define  windows  through  which  the  echo  can  be  more  easily  analyzed. 
Localized,  in  time-frequency,  noise  can  be  rejected  by  choosing  a  window  avoiding  the  noisy 
region.  For  waveforms  having  a  specific  geometric  structures  when  viewed  through  a  window, 
echo  analysis  and  noise  reduction  can  best  proceed  by  restricting  the  echo  to  this  window. 
The  formula 

PZrlxrn  =  CZLx,  x  E  CN, 

is  the  starting  point  for  studying  the  Zak  space  representation  of  the  shifts 

S%nxrn,  0  <  m  <  {R  -  1  )N. 

The  main  result  shows  that  the  restriction  of  PZRlSrnxrn  to  the  first  L  rows  of  LR  x  K 
Zak  space  is  ZLS^x,  0  <  m  <  (R  -  1  )N.  The  problem  raised  by  restricting  attention  to 
these  first  L  rows  of  LR  x  K  Zak  space  is  that  ambiguities  are  raised  by  the  inability  to 
distinguish  between  the  vectors 


S$fNxRN,  0  <r<R, 


for  fixed  0  <  n  <  N. 
For  x  E  CN,  set 


ZLx  =  [X0  X1  •  •  •  X*_ij  and  ZRLxRN  =  [K0  Yx  ■  ■  ■  YK.X) . 
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Then 
Since 
we  have 


Zrl  {Srn’xrn^  =  [DrlYk-i  Yo  ■  ■  ■  Yk- 2]  • 

CXk  —  PYk,  0<k<  K, 

PZRL  (l Srnxrn )  =  [(. PDrlP -1)  CXK—\  CXo  •  •  •  CXK_2] 


where 


PDrlP-1 


[  tr-1Dl  J 

As  a  consequence  the  restriction  of  PZrl  ( Srn~x.rn )  to  the  first  L  rows  of  RL  x  K  Zak  space 
is 

Zl  (Snx)  =  [Dl  Xk _!  X0  •••  A*_2] . 

The  relationship  between  the  remaining  rows  of 


PZrL  (S%nxrn)  and  ZL{S%x) 

is  more  complicated  due  to  the  lack  of  commutativity  between  the  FT  and  shift  operations. 
In  Section  12  we  develop  this  relationship,  but  as  yet  the  results  have  not  proved  useful  in 
echo  analysis. 

Set 

Wo  =  {(l,k):0<l<L,0<k<K}. 

For  y  €  CRN 

WqPZrlY 

is  the  L  x  K  image  formed  by  restricting  PZrlY  to  Wq.  WqPZrlY  is  a  view  of  y  through 
the  time-frequency  window  Wo- 

W0PZrl  (S%nxrn)  =  ZL  (S»  , 

and 

W0PZrL  (S%+nxrn)  =  Wo  PZrl  (S%nxrn)  ,  0  <  m  <  (R  -  1  )N. 


11  Windowed  Echos  of  Zero-padded  Discrete  Chirps 

Suppose  x  G  C'v  is  a  discrete  chirp.  In  Section  12  we  show  that  the  set 

{PZHL(5WxK7V):0<r<iV} 
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is  not  orthogonal,  but  the  set 

WQPZRt  (Vm)  =  (SrNx) ,  0  <r<N, 

is  orthogonal  and  the  expansion 

W0PZRLe  =  £  (E  a(r  +  sN ))  ZLSrNx 

T= 0  \s=0  / 

of  an  echo  e  of  xRN  is  an  orthogonal  expansion.  The  coefficients 

R—2 

^2  a(r  +  sN ),  0  <r  <  N, 

s= 0 

can  be  computed  by  inner  products.  If  the  echo  e  has  the  form 

e  =  E  a(r)5WxHiV. 

r=0 

there  is  no  aliasing  of  the  coefficients,  and  we  can  compute  the  coefficients  of  the  echo  from 
the  orthogonal  expansion 

W0PZRLe  =  f)  a(r)ZLSrNx , 

r— 0 

by  inner  products.  As  in  Section  6  dechirping  can  simplify  echo  analysis. 

An  echo  of  the  form 

e  =  £  a(r)S£f°'Vu' 

r=0 

looks  the  same  through  the  window  W0  for  all  integers  r0.  r0  can  be  determined  from 
application,  from  an  estimate  of  the  starting  time  of  the  echo  and  potentially,  from  several 
size  FZT  computations  of  the  echo. 

r(m 


Define  the  N -•periodic  multiplication  dechirping  operator  Mx.  ;  by 

0 


M^N]  =  Ir®Mx 


Mx. 


0 


Mx. 


Write 

Then 


m  =  r  +  sN ,  0  <r  <  N,  0  <  s  <  R  —  1. 


i\ARNl  -vRN  —  QsN  Qr  -vRN 

iwx*  & RN x  —  ,3iW,ux*  ^RN*- 

N  replaced  by  RN  and  L  replaced  by  RL  implies 

Zrl  (Ssrnn y)  =  DsrllZrl y, 
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and 

Recall 


W0PZRL  (Ssrnn y)  -  WoPZrly,  y  G  CRN. 


pd%lp-'  = 


tsIl 


T  =  e?mRL . 


Combining  these  results  we  have 

[M[rn)  S^NxRN)  =  W0PZRL  (M^N)SrRNxRN)  . 

Since  0  <  r  <  N, 

M^N)SrRNxRN  =  SrRN  ((S^x*)  xf , 

and 

W0PZRL  (MiRN)S%NxRN)  =  ZL  (Mx*SrNx) . 

Consider  the  echo 

{R-1)N-1 

e  =  £  a(m)S^NxRN. 

m= 0 

Dechirping  e  and  operating  by  W0PZRL  we  have 


W0PZn  (M<?">e)  =  L  £  (y  a(r  +  sW))  *(-r)*$(iW(tf). 

r=0  \s=0  / 

Special  echoes  can  be  handled  exactly  as  before  in  Section  6. 

The  Wo-window  is  emphasized  in  this  report  but  there  are  other  windows  with  potential 
applications. 

12  Linear  Shifts 

In  this  section  we  study  the  relationship  between 

PZRL  (S%nxrn)  and  ZL(S%x) 

over  complete  RL  x  K  Zak  space.  The  results  are  technical  and  have  as  yet  not  been  applied 
to  echo  analysis.  The  presentation  will  be  brief  with  several  derivations  omitted. 

In  Section  10  we  showed  that 

PZRL  (Srnxrn)  =  \DCXK-i  CX0  •  •  •  CX*_2] , 
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where  D  is  the  RL  x  RL  diagonal  matrix 

Dl 


D  =  PDrlP-1  = 


tDl  0 


0 


tr-1Dl  J 


t  =  e2mRL . 


For  the  remainder  of  this  section, 

m  =  r  +  sN,  0  <  r  <  N,  0  <  s  <  R, 


and 


r  =  t  +  uL  =  a  +  bK ,  0  <  t,  b  <  L,  0  <  u,  a  <  K. 

For  0  <b<L,  define  the  vector  ff  £  CL  by 


fL  — 


lb 

QL-b 


and  the  L  x  L  circulant  matrix  C(b)  by 

C(b)  =  F(L)D{lt)F(L)-\ 

The  difficulty  in  studying  the  shifts 

S%nxrn,  0  <m<(R-  1)N, 


by  their  RL  x  K  Zak  space  representations  is  due  to  the  noncommutativity  of  Dl  and  C. 
Dl  and  Sl  do  not  commute! 

Define  the  vector  gR  £  CL  by 


0  <  b  <  L. 


The  vectors  Mk(b)  £  CL, 

Mk(b)  =  C{b)DbLXk,0  <  k  <  K, 


are  the  blurring  factors. 

Suppose  x  £  CN  is  a  discrete  chirp.  Expanding 

C(b)  =  jEDcL(lLs  (i1)')  Dl°, 

Mk(b)  =  x(k)y/v-(b-e)tDcL  l1 

c=0 


and 
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and  placing  this  result  in 

PZRL  (SrRNxRN)  =  Bb  [DCX*_a  •  •  •  DCVi  CX0  •  •  •  CXk^I  . 

The  term 

zCML(fl,  b ) 

is  the  problem.  It  is  responsible  for  obstructing  the  construction  of  a  zero-padded  discrete 
chirp  whose  linear  shifts  form  an  orthogonal  set.  In  fact  we  can  show  that  the  inner  product 

(PZRLxRN,PZRL(SRNxRN)) 

is  not  zero.  This  inner  product  is 

(C ZLx,  CZL  (SNx))  +  (CZLx,  zCMi(l,  0)). 

The  first  inner  product  vanishes.  Since 

Mz,(l,0)=  [m*-i(1)  0l  •••  0L], 

the  second  inner  product  is 

f>c(X0,M*_a(l))  =  —R{X0, 

c= 0 

Prom 

MK-i(  1)  =  vx(l)lL, 

we  have 

(X0,MK-1(1))  =  Lv-1x*(1), 

and 

(• PZRLxRN,PZRL  (Srnxrn))  =  -RLv~lx*{  1)  ^  0. 


13  Filtering  Algorithms 

The  algorithms  of  Section  11  are  based  on  windowing  the  RL  x  K  Zak  transform  of  an 
echo  e.  In  this  section  a  more  complex  filtering  of  this  Zak  space  image  will  be  introduced. 
Surprisingly,  perhaps,  the  resulting  algorithm  has  simple  structure.  The  range  of  applications 
is  the  same  as  that  handled  by  the  Vbo-window  algorithm.  The  eventual  goal  is  to  combine 
these  approaches  in  a  multiple  chirp  interrogation  of  a  general  target  scene. 

For  y  G  CRN,  write 

y  =  [ys]0<s<R.  y^cw,  o<s<i?. 

The  basic  idea  underlying  the  approach  in  this  section  is  to  segment  an  echo 

RN- 1 

y=  E  a(m)S%NxRN, 

m=0 


x  a  discrete  chirp, 
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into  the  R  segments 


ys  e  CN,  0  <s<R, 

and  then  to  compute  the  FZT  of  each  segment 

ZLys,  0  <s<R. 

Under  certain  conditions  the  coefficients  of  y  can  be  computed  by  first  computing  the  sums 

ZLys  +  ZL  y8+1,  0  <s<R, 

and  then  applying  the  methods  of  the  preceding  sections. 

Suppose  y  e  CRN .  The  first  result  relates  Zrl y  to  the  FZT  of  the  segments  of  y 

zLys,  o  <  s  <  R. 


First  observe  that 

MRLy  =  [MLyr]0<r<R . 

and 

PZRL y  =  PF(RL)Mrl y  =  PF(RL)  [ML yr}0<r<R  •  (2) 

The  RL  CT  FFT  factorization  will  be  used  to  compute  PZRRy  in  terms  of  ZR ys,  0  < 
s  <  R. 

Set 

C*  =  ^  [hC-1  •  •  •  C'-M] 

and 

F  =  (FiR)-1  0 IL)  C*. 

F  can  be  viewed  as  a  Zak  space  filtering  operation. 

As  before  we  write 


m  =  r  +  sN,  0<r<A',  0  <  s  <  R, 


and 


r  =  a  +  bK  =  t  +  uL,  0  <  a,  u  <  K,  0  <b,t<L 
and  suppose  x  is  a  discrete  chirp  satisfying  the  L  x  K  Zak  space  condition. 

F PZRL  (S%nxrn)  =  (. F(R r1  0 IL)  ( BsLZl  (SrNx)  +  zD sLML(a,  b))  . 


Since 

and 


(F(R)-1  ®  IL)  DsL  =  ef  0  IL 
(F(R)-1  ®  IL)  (zDsL)  =  -ef  0  IL  +  ef+1  0 IL. 


we  have  the  following. 
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Suppose  y  is  an  echo  of  the  form 


y=-£a(r  +  sN)S$Nxm, 


for  fixed  0  <  s  <  R.  Then 


Then 


ZLys  =Y,a(r  +  sN )  (ZL  ( SrNx )  -  ML(a,  b )) , 


ZLys+\  =  E  a(r  +  sN)ML(a,  b ) 


ZLys>  =  o^,  s'  ^  s,  or  s'  ±  s  +  1. 


E  a(r  +  s7V)Zl  (SrNx)  =  ZLys  +  ZLys+1 


and  we  can  compute  the  coefficients  of  y  by  the  previous  methods  from  the  sum 

ZlY s  +  ZLys+1. 

The  same  argument  shows  that  for  echoes  of  the  forms 

•  y  =  ot(r  +  sN)Srj^NxRN ,  for  fixed  0  <  s  <  R. 

•  y  =  T*t~?a(r  +  s1N)SrI^lNxRN  +  '£,r=oair  +  s2N)S^2Nxrn ,  for  fixed  0  <  slt 
s2  <  R,  with  |sj  -  s2\>2 


Ri  N- 1 

y  =  EE  a  (r  +  SjN)  S'^’]vS',iVxi?Ar,  for  fixed  0  <  si  <  s2  <  ■  ■  ■  <  <  R 

j= 1  r=0 


where  |sj  —  sJ+i|  >  2,  1  <  j  <  R\, 

the  coefficients  of  the  echoes  can  be  computed  by  the  methods  of  the  preceding  sections  from 
the  sums 

ZlYs  +  ZLys+i ,  0  <  s  <  R. 

However  if  echo  y  has  the  form 

y  =  E  a{r)sRN*-RN  +  E  a(r  +  n)sr+nnxrn> 


ZLy i  =  E  <*(r)ML{a ,  b)+Y  a(r  +  N)  ( ZL  (SrNx)  -  ML(a,  b)) ,  r  =  a  +  bK , 

r=0  r=0 

from  which  we  see  that  the  information  containing  the  coefficients  of  the  first  sum  to  that  of 
the  second  sum  cannot  be  decoupled. 
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14  Material  Identification  in  Zak  Space 

In  a  previous  report  [1]  we  developed  algorithms  for  numerically  deconvolving  Krueger’s 
formula 

ea(t)=  pa(t  -  u)xTa{u)du,  teft, 

J  —  oo 

relating  the  reflectivity  kernel  pa  of  a  dielectric  material  to  the  echo  ea  resulting  from  the 
interaction  of  an  orthogonally  incident  chirp  x ^  with  the  material.  The  goal  in  this  section 
is  to  describe  a  new  deconvolution  algorithm  using  the  Zak  space  framework  developed  in 
the  preceding  sections. 

Assume  throughout  that  pa  is  supported  in  an  interval  [0,  J),  J  estimated  in  some  range 
but  not  known  and  x %  a  waveform  supported  on  the  interval  [0,  T).  The  echo  ea  is  then 
supported  in  the  interval  [0,  T  +  J).  Until  required  we  do  not  assume  that  x Ta  is  a  chirp. 
Choose  an  integer  R  satisfying  the  condition 

J  <  RT. 


If  there  are  several  possible  materials  under  investigation,  we  can  choose  R  sufficiently  large 
so  that  this  condition  holds  for  all  the  estimated  J’s. 

The  expression  for  the  echo  will  be  discretized.  Choose  a  positive  integer  A"  and  set 
M  =  NR.  Approximate  the  integral  by  the  Riemann  sum  corresponding  to  the  points 

^ RT  =  ^T,  0  <m<M, 

M  N 

with  the  resulting  approximation 

Sample  the  echo  at  the  points 

~RT,  0  <  n  <  {R  +  1  )N, 

M  ~ 

with  the  resulting  approximation 


Setting 


x=  ( 


mr 

AC 


e  = 


0<m<iV 


eal^T 


N  J\o<n<(R+l)N 


and 
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the  approximation  can  be  written 


rp  M—  1 

e=^E  P(m)5(R+l)JVx 


RN 


m= 0 


e  is  an  echo  of  the  zero-padded  vector  x^+h'v,  x  G  CA' . 

In  applications  it  sometimes  useful  to  modify  the  discretization  by  sampling  the  integral 
at  the  points 


m 

N 


T,  1  <  m  <  M, 


and  the  echo  at  the  points 


n 

N 


T,  1  <n<(R+  1)N, 


with  the  resulting  approximation 

rp  M 


«■  -  n  S '  (Sr)  (" irT )  •  i  s » s  (« + ^ 

By  the  change  of  variables 


ea 


n'  +  l, 
N 


m'  =  m  —  1  and  n'  =  n  —  1, 

T^1  /m'  +  l  \  rrfn'-m  ■ 


T  =  77  £  a. 


/V  ^ 

m'=0 


N  ~TiXa 


Setting 


and 


we  have 


[  T  frnm\ 

r  (  T\ i 

K  UT)J 

,  e  = 

e“  VN  ). 

0  <m<N 

N 


\<n<{R+l)N 


p{m)=  Pa(jjT^j 


Km<M 


rp  M  —  1 

e  =  T7  £  +  1)S$+1)nx.(r+1)n. 


m—0 


The  second  sampling  is  preferred  whenever  we  want  to  avoid  pa  at  the  origin. 

Methods  for  analyzing  such  echoes  in  Zak  space  were  discussed  in  the  preceding  section 
when  x J  is  a  chirp 


= 


e^  e2^*,  0  <  t  <  T, 
0,  otherwise, 


t  G 


such  that 


N  =  7T2 

is  a  positive  integer  and  the  periodic  modulo  N  condition  holds. 


emNe2mPt  _  ^ 


Assume  for  the  rest  of  this  report  that  a  chirp  satisfies  these  two  conditions. 
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Suppose  now  that  xj  is  a  chirp  in  the  formula 

rp  (Rl-l)JV-l 

e=U  E  Ki  =  «+1, 

■'V  m= 0 

derived  above.  Set  N  =  LK.  Windowing  to  the  first  L  rows  of  R\L  x  K  Zak  space,  we  have 
that  the  expansion 

rp  {Ri—1)N—1 

W0PZRlLe  =  -  pMZlS%x,  P  =  P(R1L,  R1), 

^  m=0 

is  an  orthogonal  expansion  after  collecting  like  powers  of  Sn  and  we  can  compute  the  sums 

iii-i 

E  p(r  +  sA/'),  0  <  r  <  N, 

8=  0 

by  inner  products. 

For  sufficiently  large  T  we  can  in  the  absence  of  noise  compute  the  scattering  coefficients 
themselves.  Suppose 

J  <T 

and  we  take  R=  1.  Windowing  to  the  first  L  rows  in  2 L  x  K  Zak  space 

W0P(2L,  2 )Z2Le  =  ^  E  pMZlS^x 

771= 0 

is  an  orthogonal  expansion  and  we  can  compute  the  p(m),  0  <  m  <  N,  by  inner  products. 

Figures  17  -  19  illustrates  this  orthogonality.  In  Figure  19,  usual  matched  filter  processing 
result  is  also  displayed  and  compared.  Note  that  the  matched  filter  processing  is  over  the 
entire  space  and  orthogonality  is  lost.  The  relevant  parameters  for  these  figures  are  J  =  10, 
T  =  20,  R  =  1,  7  =  .25  and  /  =  20.  The  Zak  transform  parameters  are  L  =  10  and  K  =  10. 
Figures  20  and  21  display  equivalent  information  in  the  presence  of  noise. 

Figure  17:  Test  reflectivity  kernel 


Suppose  N  has  factorization  N  =  LK  such  that 


J  < 


K' 
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Figure  18:  Sampled  echo  from  discrete  chirp  pulse  upon  interrogating  the  test  material 


real  part  imaginary  part 


Figure  19:  Recovered  reflectivity  kernels 


windowed  inner  products  matched  filter  processing 


In  2L  x  K  Zak  space 

W0P(2L,2)Z2Le  =^E  P(m)ZLS%x. 

771=0 

The  supports  of  the  terms  of  the  summation 

p(m)ZLS^x,  0  <  m  <  L, 

are  pairwise  disjoint  lines  in  Lx  K  space.  |p(m)|  is  the  common  absolute  value  at  each  point 
of  the  line 

Ex(m),  0  <  m  <  L. 

For  a  noisy  echo  we  can  estimate  \R(rn)\  from  the  possibly  varying  absolute  values  on  the 
line. 

For  the  chirp  pulse  the  integral  is  approximated  by  a  Riemann  sum  corresponding  to  the 
points 

-tjRT  =  -~T  =  m~,  0  <m<M. 

M  N  7  T 

The  larger  we  take  7 T,  the  finer  the  sampling  grid,  and  the  better  the  approximation.  The 
same  is  true  for  the  reflectivity  kernel  samples  p(m),  0  <  m  <  M. 

Figures  23  and  24  are  included  to  illustrate  how  the  Zak  space  representations  reflect  the 
effective  supports  of  the  reflectivity  kernels. 


•  FA9550-05-C-0029:  Final  Executive  Summary 


38 


Figure  20:  Sampled  noisy  echo 


Figure  21:  Recovered  reflectivity  kernel  using  windowed  inner  products 


Simplest  denoising  in  Zak  domain  is  used  to  estimate  the  reflectivity  kernel  coefficients 
in  the  presence  of  noise.  The  number  of  lines  in  the  Zak  space  representation  of  the  echo  is 
determined,  then  the  values  along  each  line  are  averaged.  Figures  25  -  27  display  the  results 
of  denoising.  Zak  transform  parameters  are  L  —  40  and  A  =  40. 

Setting  L  =  20,  we  have  another  estimate  for  the  first  20  coefficients  of  the  reflectivity 
kernel.  The  comparison  between  these  estimates  are  displayed  against  the  test  kernel  in 
figure  32.  The  first  plot  is  that  of  the  estimate  from  using  L  —  40.  The  second  plot  is  the 
estimate  from  using  L  =  20.  The  last  is  the  average  of  the  two  sets  of  the  estimates. 

There  is  a  second  approach  which  potentially  increases  the  accuracy  of  the  integral  ap¬ 
proximation  and  provide  better  resolution  for  the  reflectivity  kernel  samples.  As  before 

J  <RT. 


Choose  any  positive  integer  A  and  set  M  =  RNK.  Sampling  the  integral  and  the  echo  at 
the  points 


—  RT  —  n  T 
M  KN  ’ 


we  have  the  approximation 


T  (  m  \  T  (n  —  m  \ 

AW  Pa  \KN)  Xa  l  KN  T)  ’ 


0  <  n  <  (R  +  1)  AW 
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Figure  22:  Recovered  reflectivity  kernel  using  matched  filter  processing 


real  part  imaginary  part 


Figure  23:  Test  reflectivity  kernel  and  the  Zak  transform  of  the  echo 


real  part  imaginary  part  absolute  value 


The  sequence 


where  x „  is  a  chirp  is  the  oversampled  chirp  in  Section  2.  Future  work  will  study  the  Zak 
space  structure  of  oversampled  chirps,  shifts  and  echoes. 

15  Atmospheric  Effects 

15.1  Introduction 

In  this  section  we  begin  a  preliminary  investigation  into  the  effects  of  atmospheric  interference 
on  a  waveform  interrogating  a  ground  based  dielectric  material  from  a  space-based  radar.  We 
propose  two  multiple  measurement  strategies  for  removing  this  interference  by  singular  value 
decomposition  methods.  In  future  work  we  plan  to  study  methods  based  more  significantly 
on  the  electromagnetic  differences  between  the  atmosphere  as  a  magnetically  biased  plasma 
and  an  interrogated  dielectric  material. 
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Figure  24:  Test  reflectivity  kernel  and  the  Zak  transform  of  the  echo 


real  part  imaginary  part  absolute  value 


Figure  25:  Noise-free  echo 


real  part  imaginary  part 


Throughout  this  section  atmospheric  interference  is  modeled  as  a  linear  convolution. 
Suppose  a  pulse  xa  is  transmitted  from  a  space-based  radar,  passes  through  atmospheric 
interference  modeled  as  a  linear  convolution  with  an  impulse  response  function  hi,  inter¬ 
rogates  a  dielectric  material  having  reflectivity  kernel  pa  and  returns  through  atmospheric 
interference  now  modeled  by  an  impulse  response  function  hi-  Suppressing  signal  translation 
due  to  travel  time  the  echo  ea  is  the  triple  convolution 

h2(w)  J  pa{v ) 

By  Fubini’s  theorem  we  can  under  natural  assumptions  on  the  waveforms  in  the  triple  con¬ 
volution  interchange  orders  of  integration  and  write 


J  h\{u)xa{t  —  w  —  v  —  u)dudvdw. 


ea(t )  =  J(h2*  hi){r){pa  *  xa)(t  -  r)dr, 


where 


h(r)  —  ( h2  *  hi)(r)  =  /  h2(w)hi(r  —  w)dw 
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Figure  26:  Noise-added  echo 


real  part 


imaginary  part 


Figure  27:  Test  and  estimated  reflectivity  kernel 


and 


Pa  *  xa(t)  =  J  pa(v)xa(t  -  v)dv. 


Suppose 


supppa  C  [0,  T) 

and  that  R\  and  R2  are  integers  such  that 

supppa  C  [0,  RiT) 


and 


supph  C  [0,  R2T). 

The  echo  ea 

is  supported  in  the  interval 

[0,  RT ),  R  =  Ri  +  i?2  +  1 

and 

rR'iT 

ea(t)  =  /  h(r)  (pa  *  xa)  ( t  -  r)dr, 

Jo 

t  €  ft, 

where 

fRiT 

(Pa  *  Xa)  ( t )  =  /  pa(u)xa(t  -  u)du , 

Jo 

t  e  ft 
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Figure  28:  Noise-added  echo 


real  part  imaginary  part 


Figure  29:  Test  and  estimated  reflectivity  kernel 


real  part 


imaginary  part 


Choose  a  positive  integer  N  and  set 

Mi  =  RiN,  M2  =  R2N  and  M  =  RN. 
Discretizing  the  integral  formula  for  the  echo  by 


and 


m  „  nr 

r  =  W2R2T=NT’ 


t  =  W2R2T  =  WT’  0£n<M- 


we  have  the  approximation 

rp  M2  — 1 


>(W  »*»<*• 


Discretizing  the  integral  formula  for  pa  *  xa  by 


u  =  lkR2T=TiT ’ 


we  have  the  approximation 


n  —  {m  +  k) , 

7T 


T  ,  0  <  n  <  RN, 


*t 
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Figure  30:  Noise-added  echo 


real  part  imaginary  part 


Figure  31:  Test  and  estimated  reflectivity  kernels 


Using  the  notation  of  the  preceding  sections  with 


and 


r  /  n  _M 

e  = 

e“  (ivT). 

,  x  = 

0 <n<RN 

- 1 

h 

B 

C3 

p(k)  =  Pa  (JyTj  ,  0  <  k  <  Ml,  h(m)  =  h  (£r)  ,  0  <  m  <  M2 


we  have 

rp  2  M2  —  1  Mi  — 1 

e  =  M2  £  tl  h{m)p{k)S^+kxM. 

iV  m= 0  fc=0 

The  echo  e  is  a  linear  combination  of  shifts  of  xM ,  however  the  coefficients  are  corrupted  by 
the  samples  of  h. 

In  order  to  place  these  formulas  Within  the  SVD  framework  we  define  the  sequences 


e(n)  =  ea  ,  n  G  Z,  x{n)  =  xa  ,  n  G  Z, 

P(*)  =  Pa  ,  k  G  Z,  /r(m)  =  h  (%t)  ,  m  6  Z. 


and 
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Figure  32:  Test  and  estimated  reflectivity  kernels 


Each  sequence  is  causal  and  has  finite  support.  The  echo  e  is  the  triple  linear  convolution 
of  sequences 

T2 

e  =  ]y2  h*(P*x)- 

The  problem  is  to  estimate  p  *  x  from  e.  Two  multiple  measurement  strategies  will  be 
outlined  for  removing  h  by  SVD  methods.  These  methods  are  described  in  detail  in  the 
following  sections  and  codes  have  been  implemented.  The  extensive  treatment  of  these 
SVD  methods  can  potentially  lead  to  other  multiple  measurement  strategies  and  multiple 
waveform  strategies  for  removing  atmospheric  effects.  These  strategies  will  be  studied  in 
future  work. 


15.2  Multiple  Measurement  Strategies 

15.2.1  Strategy  1 


We  propose  two  transmissions  of  the  same  pulse  xa  from  a  space-based  radar  with  the 
assumptions  that  the  atmospheric  interference  h  is  the  same  for  both  transmissions,  but  the 
ground  region  under  interrogation  is  a  known  dielectric  material  pi  for  the  first  transmission 
and  an  unknown  dielectric  material  p2  for  the  second  transmission.  The  corresponding  echoes 
are 

e\  —  h  *  (pi  *  x)  and  e2  =  h  *  (p2  *  x). 


Set 


Si  =  pi  *  x  and  s2  =  p2*  x. 


The  problem  is  to  compute  s2  with  si  known. 


15.2.2  Strategy  2 

We  propose  two  transmissions  of  the  same  pulse  xa  from  a  space-based  radar  interrogating  the 
same  dielectric  material  p ,  but  the  first  passing  through  atmospheric  interference  h  =  hi  *  hi 
and  the  second  passing  through  atmospheric  interference  h!  =  [h\  +  A)  *  (h\  +  A) ,  where  A 
is  a  known  perturbation.  The  corresponding  echoes  are 


e  =  h  *  (p  *  x),  e'  —  h!  *{p*  x). 
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The  problem  is  to  estimate  h.  As  in  Strategy  1  if  h  is  computed,  we  can  then  compute  p*x. 

There  are  several  potentially  useful  variants  of  measurement  strategies  1  and  2.  Different 
pulses  can  be  transmitted,  pulse-trains  suitably  separated  in  time  can  be  transmitted.  Angle 
of  incidence  dependence  can  be  used  to  distinguish  between  h  and  p. 


16  Future  Efforts 

To  realize  the  potential  of  the  methods  and  algorithms  developed  in  this  report  as  tools 
for  material  identification  we  will  carry  out  research  leading  to  algorithms  in  the  following 
topics. 

•  Refinement  and  extension  of  existing  algorithms. 

•  Waveform  design  in  Zak  space. 

•  Multi-measurements  of  multiple  targets,  including  the  dependence  of  reflectivity  kernels 
on  incident  angle. 

•  Atmospheric  effects. 

16.1  Refinement  and  extension  of  existing  algorithms 

We  identify  the  following  topics  for  future  work 

•  Varying  size  FZT  representations. 

•  Windows  in  Zak  space. 

•  Oversampled  chirp  pulses. 

Throughout  this  report,  a  fixed  factorization  N  =  LK  is  taken  to  represent  a  discrete 
waveform  (before  zero-padding)  of  size  N  in  L  x  K  Zak  space.  For  each  factorization  a 
noise-free  echo  of  a  discrete  chirp  in  Zak  space  is  supported  by  a  series  of  parallel  lines  of 
slope  1,  but  the  aliasing  depends  on  the  factorization.  We  will  study  how  processing  the 
echo  relative  to  several  factorizations  can  be  used  to  remove  this  aliasing. 

In  numerical  experiments  we  have  seen  that  echo  reconstruction  from  noisy  echoes  varies 
as  the  factorization  varies.  It  seems  that  noise  is  distributed  differently  for  different  factor¬ 
izations  of  N  that  is  in  varying  size  Zak  space  representations.  We  will  study  this  dependence 
to  improve  on  noise  reduction  algorithms. 

In  this  report  a  specific  Zak  space  window  Wo  was  specified  so  that  through  Wo  the  col¬ 
lection  of  linear  shifts  of  a  critically  sampled  zero-padded  chirp  pulse  is  orthogonal.  However 
information  is  lost  through  the  window.  The  resulting  windowed  orthogonal  expansion  of 
the  echo  is  aliased.  We  plan  to  study  new  windows  either  subsets  of  Zak  space  or  unitary 
transforms  of  windows  in  Zak  space  for  removing  this  aliasing  for  critically  sampled  chirp 
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pulses  and  more  general  sampled  waveforms.  These  new  windows  can  also  be  matched  to 
the  noise  characteristics  in  a  specific  application. 

Oversampled  chirp  pulses  were  mentioned  in  the  report  with  a  formula  relating  oversam¬ 
pled  chirp  pulses  to  critically  sampled  chirp  pulses.  We  plan  to  study  the  extension  of  the 
Zak  space  framework  to  oversampled  chirp  pulses.  This  effort  will  provide  greater  flexibility 
in  the  use  of  Zak  space  methods  to  multiple  chirp  pulses  of  varying  time  duration  and  chirp 
rates. 


16.2  Waveform  design 

Zak  space  methods  are  specially  tuned  to  critically  sampled  echoes  of  chirp  pulses  providing 
a  windowed  representation  supported  on  series  of  parallel  lines.  We  plan  to  study 

•  The  representation  of  echoes  resulting  from  more  general  waveforms,  including  multiple 
critically  sampled  chirp  pulses,  upchirps  and  downchirps. 

•  Waveform  design  directly  in  Zak  space  including  lines  at  varying  slopes,  polygonal  and 
other  geometric  structures. 

In  the  second  study  the  inverse  Zak  transform  will  be  used  to  construct  the  appropriate 
discrete  waveform  and  several  methods  will  be  used  to  construct  continuous  waveforms  and 
sampling  rates  supporting  these  one-dimensional  discrete  waveforms. 

We  are  especially  interested  in  designing  Zak  space  waveforms  which  behave  well  with 
respect  to  shifts  (lines)  and  rotations  (accelerated  targets).  Since  the  Zak  space  representa¬ 
tion  of  FT  is  essentially  90°  rotation,  geometric  structures  invariant  under  90°  rotation  may 
be  of  special  interest.  In  a  previous  report  [1]  we  use  this  idea  to  construct  an  orthogonal 
eigen  vector  basis  of  the  FT  having  minimal  time-frequency  support. 


16.3  Multiple  measurements/multiple  targets 


In  [1]  we  developed  and  tested  algorithms  for  distinguishing  I0  stationary  dielectric  targets 
along  with  their  reflectivity  kernels  using  Jo  strip  SAR  measurements  from  chirp  pulse  echoes 
under  the  assumption  that  the  reflectivity  kernels  are  independent  of  the  incident  angle  of 
the  interrogating  chirp  pulses.  The  results  of  Section  14  will  be  used  to  develop  a  similar 
algorithm  in  the  Zak  space  framework. 

As  in  Section  14  suppose  x Ta  is  a  chirp  pulse 


x 


T 

a 


eiri-yt2  e2iri/3t ,  g  <  t  <  T,  f  £  p 

0,  otherwise  "  ’ 


satisfying  N  =  7 T2  is  a  positive  integer  and 


(FiNe2mf  =  1,  f  =  pT. 


Assume  throughout  that  the  relaxation  constant  J  of  all  the  dielectric  materials  interrogated 
satisfy 


J  <  T. 
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The  formula  for  the  windowed  critically  sampled  echo  in  Section  14  does  not  include  the 
shift  due  to  the  time  traveled  between  antenna  and  target.  Suppose  a  chirp  xj  having  carrier 
frequency  f3  is  transmitted  from  antenna  at  v  and  orthogonally  interrogates  a  dielectric 
material  with  reflectivity  kernel  pa  at  u.  The  echo  at  v  is 


fT 

5b(*)  =  /  Pa(t)xa(-U0  +  t  -  v)dv, 
JO 


where  u0  =  -r\v  —  u|.  Write 


.  j- 

U0  =  iQ-jy  +  CO, 


where  l0  >  0  is  an  integer  and  e  <  jj.  Then 

T 


a  (t  +  lo—'j  =  J  pa(v)xl(-eo  +  (t  -  v))dv. 


Sampling  the  integral  at 


777  77, 

v  =  —  T,  0  <  m  <  N,  and  t  =  —T,  0  <  n  <  2 N, 


we  have  the  approximation 
'n  +  hr 


eo 


N 


rKI 


If  eo  =  0,  we  can  argue  as  before  to  compute  the  reflectivity  kernel  samples.  We  assume  the 
model  introduced  in  2  that 

°  <  e0  < 

and  that  for  e0  >  0,  we  can  make  the  approximations 

xTa(-to  +  T)  =  xTa(0)e~^ 


and 


£  (~eo  +  =  xl  (£t)  1  <n<N. 


Setting 


e  = 


e  ■  n  +  l°T 
al  TV 


J  0<n<2n 


>  p(m)  =  P  (jjT)  >  0  <  m  <  N, 


and 


x(e0)  = 


zli-eo  +  T) 
(_eo  +  jjT) 


=  uT2wi/3£0x,  x  =  x(0), 


$  •  FA9550-05-C-0029:  Final  Executive  Summary 

) 

we  have 


48 


e  =  pHS?/1  (Spx) 


2  N 


m= 0 


and 


W0P(2L,  2 )Z2Le  =  £  p(m)ZLS%x. 


m=0 


Suppose  a  chirp  xj  having  carrier  frequency  j3j  is  transmitted  from  an  antenna  at  Vj  and 
interrogates  dielectric  materials  having  incident  angle  independent  reflectivity  kernels  /$)  at 
Uj,  0  <  i  <  Io-  Set 

2  T 

ui,j  =  ~  I  vj  I  =  jy  ^ 


c 

T 


where  Iq  >  0  is  an  integer  and  Cij  <  jj,  0  <  i  <  Io- 
Setting 

'n  +  l0r 


ei 


A’ 


,  x,-  = 


J  0<n<2W 


0  <n<N 


and 


PWM  =  p«  (~T)  ,  0  <  m  <  N, 

we  have  from  the  preceding  discussion  that 

M/„P(2L,  2)Z2Ie2  =  A  £  (z  2iS^x2. 


m— 0  \  i=0 


By  orthogonality  we  compute 


/o-i 


Aj  (m)  =  Yl  P{i)  ,  0  <m  <  N. 

i= 0 


In  matrix  notation 


[Aj(m)]0<m<jy  —  REj, 


where  i?  is  the  N  x  I0  matrix 


and  Ej  is  the  vector  C/o. 


R  =  [pw(m)] 


0 <m<N,  0<i<lo  ’ 


Ej  =  [c-2^^] 

Processing  over  Jo  antenna  positions  we  compute 

A  =  RE, 


0<i<Io 


where 


A  —  [^j(m)]o<m<W,0<j<^o  ’ 
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and 

E  =  [Eq  Ei  ■  ■  ■  EJo-i\ . 

E  models  the  geometry  of  the  measurements.  In  [1]  we  estimated  R  by  pseudo-inversion  of  E 
and  tested  the  accuracy  of  this  inversion  for  several  model  geometry  matrices  E  depending  on 
the  number  of  antenna  positions  J0  the  number  of  dielectric  targets  J0,  the  spacing  between 
antenna  positions  and  the  spacing  between  target  positions. 

Suppose  that  the  reflectivity  kernels  depend  on  the  incident  angle.  The  Nxlo  reflectivity 
kernel  matrix  R  in  the  above  discussion  is  replaced  by  Jo  N  x  /0  reflectivity  kernel  matrices 

R(j),  0  <j<J0, 

and  we  can  compute  by  the  above  processing  the  J0  vectors  in  CN 

[Ai(  rn)]o<m<N  ~  R(j)Eji  0  <  j 1  <  Jb. 

Ej  is  a  vector  in  C/o  depending  on  the  carrier  frequency  of  Xj  and  on  the  distances  between 
antenna  position  Vj  and  Iq  dielectric  targets. 

We  plan  to  develop  algorithms  for  computing  the  matrix  Rj  by  replacing  the  chirp  pulse 
Xj  by  a  sequence  of  chirp  pulses  at  different  carrier  frequencies  sufficiently  separated  in  time 
such  that  the  echo  for  each  of  the  chirp  pulses  in  the  sequence  can  be  independently  measured 
and  processed.  This  has  the  effect  of  replacing  the  vector  Ej  in  R(j)Ej  by  a  Iq  x  S'  matrix, 
S  the  number  of  chirp  pulses  in  the  sequence.  The  key  to  this  approach  is  to  choose  the 
different  carrier  frequencies  so  that  the  pseudo-inversion  of  the  resulting  I0  x  S  matrix  will 
lead  to  a  good  estimate  of  R(j). 

The  computation  of  R(j)  provides  an  estimate  of  the  reflectivity  kernels  of  the  I0  dielectric 
targets  at  incident  angles  determined  by  the  geometric  relationship  between  antenna  position 
Vj  and  the  7o  dielectric  targets.  Processing  over  several  antenna  positions  we  get  an  incident 
angle  dependent  collection  of  reflectivity  kernel  estimates  which  we  can  compare  with  existing 
libraries. 

In[l]  we  propose  several  models  of  incident  angle  dependence  for  the  purpose  of  study¬ 
ing  the  effects  of  these  models  on  algorithms  for  distinguishing  a  multiple  dielectric  target 
environment  from  multiple  SAR  measurements.  For  example  if  an  antenna  at 

v  =  (0  ,v,h) 

transmits  a  pulse  chirp  which  interrogates  a  dielectric  target  at 

u  =  (x,y,  0) 

we  assume  that  the  dependence  on  the  reflectivity  kernel  at  u  on  incident  angle  could  be 
expressed  by 

(x,0,-h) 

’  —h)\ 

In  words  the  reflectivity  kernel  depends  upon  the  slant  direction  and  not  on  the  broadside 
direction  between  target  and  antenna.  The  justification  for  this  model  is  that  the  antenna 
footprint  is  relatively  narrow  in  the  along  track  direction  so  broadside  variation  is  small. 
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The  key  to  future  work  using  this  approach  is  to  determine  more  realistic  models  of 
incident  angle  dependence  and  then  to  embed  these  models  in  algorithms. 

In  the  preceding  discussion  we  have  not  taken  into  account  the  possibility,  even  proba¬ 
bility,  that  several  dielectric  targets  are  equidistant  to  given  antenna  position.  If  this  is  the 
case  the  algorithms  described  above  compute  the  sum  of  reflectivity  kernels  of  the  equidistant 
targets.  This  is  a  standard  SAR  problem.  An  important  part  of  a  SAR  algorithm  package  is 
usually  devoted  to  distinguishing  equidistant  targets  from  one  antenna  position  by  utilizing 
measurements  from  other  antenna  positions.  Although  we  have  not  included  a  discussion  of 
this  method,  they  will  be  an  important  part  of  the  algorithmic  development  in  this  effort  in 
SAR  system  algorithmic  package. 

16.4  Atmospheric  effects 

In  this  report  we  study  the  problem  of  removing  the  ionospheric  effects  from  waveforms 
interrogating  ground  based  dielectric  material  from  a  space-based  radar  by  proposing  two 
measurement  strategies  which  in  principle  provide  sufficient  information  to  filter  out  the 
ionospheric  effects  by  SVD  methods.  Unfortunately  these  methods  can  be  highly  fragile  in 
noise. 

We  plan  to  study  methods  based  more  significantly  on  the  electromagnetic  properties 
of  the  atmosphere  as  a  magnetically  biased  plasma  and  especially  on  the  electromagnetic 
differences  between  this  plasma  and  the  interrogated  dielectric  material.  Either  combined 
with  SVD  framework  or  by  themselves  these  differences  will  be  the  basis  for  constructing 
filters  to  remove  ionospheric  effects.  We  intend  to  study 

•  The  Faraday  rotation  introduced  by  waveform  propagating  through  a  magnetically 
biased  plasma  (dielectric  tensor). 

•  The  model  of  magnetically  biased  plasma  as  a  lossy  dielectric  whose  electric  suscepti¬ 
bility  is  negative  [4]  in  distinction  from  a  real  dielectric  material. 

•  The  impulse  response  function  for  deterministic  ionosphere  as  presented  in  [5]. 

Dr.  Richard  Albanese  of  Brooks  City  Base  has  shown  that  the  echo  resulting  from  wave 
propagation  through  a  magnetically  biased  plasma  followed  by  dielectric  interrogation  and 
returning  through  the  plasma  has  a  measurable  Faraday  rotation  which  (in  most  cases)  is 
due  solely  to  the  plasma.  Measurements  of  these  Faraday  rotations  at  several  frequencies 
can  be  used  to  extract  the  electromagnetic  properties  of  the  plasma  at  transmission  time. 
We  intend  to  use  this  result  to  build  an  algorithm  removing  the  ionospheric  effects. 

The  deterministic  ionospheric  impulse  response  function  depends  upon  electron  density 
profiles.  However  so  far  these  profiles  cannot  always  be  relied  upon.  We  are  most  interested 
in  the  character  of  this  impulse  response  function  and  not  on  its  detailed  expression.  For 
example,  we  intend  to  use  the  fact  that  the  time  duration  of  this  impulse  response  function 
is  usually  larger  than  the  relaxation  time  of  the  dielectric  material  as  a  basis  for  removing 
deterministic  ionospheric  effects. 
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